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Abstract 

A collisionless plasma is modelled by the Vlasov-Poisson system in three space dimensions. A fixed 
background of positive charge, which is independent of time and space, is assumed. The situation in 
which mobile negative ions balance the positive charge as \x\ — > oo is considered. Hence the total positive 
charge, total negative charge, and total energy are infinite. Smooth solutions with appropriate asymptotic 
behavior for large \x\, which were previously shown to exist locally in time, are continued globally for 
spherically symmetric data. This is done by showing that the charge density decays at least as fast as 
|:r| -4 . Finally, an increased decay rate of \x\~ 6 is shown in the general case without the assumption of 
spherical symmetry. 



INTRODUCTION 



Let F : R 3 -J 
/ : [0, oo) x K 3 x 



[0,oo), /o : R 3 x M 3 -> [0,oo), and A : [0,oo) x R 3 -> M 3 be given. We seek a solution, 
R 3 — > [0,oo) satisfying 

d t f + v-V x f-{E + A)-V v f = 0, 



p(t,x) = J(F(v)-f(t,x,v))dv, 
E(t,x) = J p (t,y) J ^dy, 
f(0,x,v) = fo(x,v). 



(1) 



Here F describes a number density of positive ions which form a fixed background, and / denotes the density 
of mobile negative ions in phase space. Notice that if fo{%, v) = F(v) and A = 0, then f(t, x, v) = F(v) is a 
steady solution. Thus, we seek solutions for which f(t,x,v) —¥ F(v) as \x\ — > 00. Precise conditions which 
ensure local existence were given in |17j . It is important to notice that ([I} is a representative problem, and 
that problems concerning multiple species of ions can be treated in a similar manner. 



The paper will be divided into two main results. The first will be devoted to showing global existence 
of a smooth solution to ((T|) in the case of spherically symmetric data. This serves to continue the local 
existence result of [17] ■ The second, then, is devoted to achieving an increased rate of spatial decay of the 
charge density without the spherical symmetry assumption. We will assume throughout that F has compact 
velocity support, so that decay in v of the background is not an issue. Thus, the main difficulty in showing 
global existence arises in showing p decays rapidly enough in |x|. To see this, consider the following heuristic 



1 



argument, discussed in [T7]. Let r — \x\. Then, we typically expect the electric field, E, to decay like r 2 
for large r. If we let g = F — /, then 

d t g + v ■ V x g ~{E + A)- V v g = —{E + A) ■ V V F. (2) 

Viewing (E + A) ■ V„-F as a source term for g, we can only conclude r~ 2 decay for both g and p. But, if 
p really decayed like r -2 , the integral for E could certainly fail to decay like r -2 . In fact, we cannot show 
that g decays faster than r -2 , but due to cancellation in the v integral, p must decay faster. 

The Vlasov-Poisson system has been studied extensively in the case when F(v) = and solutions decay 
to zero as \x\ — > oo. Smooth solutions were shown to exist globally in time in [13) and independently in [11) . 
Important results prior to global existence appear in [T], [B], [7], and [TU]. Also, the method used by [13] is 
refined in [S] and [IS] . Global existence for the Vlasov-Poisson system in two dimensions was established in 
[T2] and [T5]. A complete discussion of the literature concerning Vlasov-Poisson may be found in 0]. We 
also mention [2] since the problem treated in said paper is periodic in space, and thus the solution does not 
decay for large |a;|. The works cited above make extensive use of the laws of conservation of charge and 
energy. However, in the problem considered here (and those of [3], [9], [16] . and [17]). the charge and energy 
are infinite, and it is less clear how to use the conservation laws. Therefore, the use of conservation laws 
will be an important issue, and we will utilize a lemma (stated here as Lemma 2) from [16] to deal with it 
properly. 



Preliminaries 

We will use notation which follows [17] . For q > 7 + V33, let 

, 8 
p = 4 

q 

and denote 

R(x) = R(\x\) = (l + \x\ 2 )z. 

We will use the norms 

IMloo = sup \g(z)\ 

zGR" 

\\p\\ p =\\pRV(x)\U 
\\g\\ q = \\g(l + \x\ 2 + \vn\\ L ^ Re) 

and 

lllslll = llsll, + l|Vff||,+ 

but never use LP or L q for finite p and q. We will write, for example, ||<?(i)||g for the || • || g norm of 
(x, v) i-> g(t, x, v). Notice that we may take q arbitrarily large since we take the initial data to be of compact 
support. 

Following [T2] and [T7j, we assume the following conditions hold for some C > and all t > 0, x € K 3 , 
and v G K 3 , unless otherwise stated : 




(I) F(v) = Fr(\v\) is nonnegative and C 2 with, 



(4) 



F R (0) < 0. (3) 

In addition, there is W £ (0, oo) such that 

F' R (u)<0 foru€(0,W) 
F R (u) = for u > W. 

(II) fo is C 1 with compact v-support, nonnegative, and satisfies the condition of spherical symmetry, 

fo(x,v) = f R (\x\,\v\,x-v) (5) 

which is equivalent to 

fo(x,v) = f (Ux,Uv) 

for every rotation U. 

Also, there is N > such that for \x\ > N, we have 

f (x,v) = F(v). 

(Ill) A is C 1 and 

\A(t,x)\ + \d x A{t,x)\ <CR- 2 (x) 

\V X -A(t,x)\ <CR- 4 (x). 
Furthermore, A satisfies the condition of spherical symmetry, 

A{t,x) = a(t,\x\)j^ (6) 

which is equivalent to 

A(t,x) = U T A(t,Ux) 

for every rotation U. 

Finally, we assume that there is a continuous function a: [0,T] — > R such that for |a:| > 



a(t,\x\)-^ 
\x\ 



< R 2 - p (x). 



It should be noted that the assumptions ([5]) and ([6]), imply the spherical symmetry of f(t, x, v), g(t, x, v), 
and E(t,x) for all t £ [0, oo), x £ M 3 , and v £ R 3 . Thus, where necessary we will write g(t,x,v) = 
g(t, \x\, \v\,x ■ v) when using the spherical symmetry of g. Furthermore, the spherical symmetry of E and / 
will be instrumental in the first results of the paper, while the last theorem is dedicated to eliminating the 
symmetry assumptions to conclude similar decay rates. We wish to prove the following : 

Theorem 1 Assuming conditions (I), (II), and (III) hold, there exists f £ C 1 ([0, oo) x R 3 x R 3 ) that 
satisfies (QP with \\\{F — f)(t)\\\ bounded on t £ [0, T], for every T > 0. Moreover, f is unique. 



In [17] , both local existence and a criteria for continuation of a bounded, unique solution of (JT|) are shown. 
Thus, to prove Theorem 1 we will need to establish the continuation criteria for all T > 0. Specifically, 
Theorems 2 and 3 of [T7], when combined, state : 

Theorem 2 Assume q > 7 + \/33 and conditions (I), (II) and (III) hold, without $5$ and (0|). Let f be a 
C 1 solution of (Qp on [0,T] x I 3 x I 3 with T > 0. If || j(F - f)(t) dv\\ p is bounded on [0,T], then we may 
uniquely extend the solution to [0, T + 8] for some 6 > with \\ \(F — f)(t)\ \\ bounded on [0, T + S\. 



Therefore, to prove Theorem 1, we will find a solution to (jTJ) on [0,T] for some T > using the local 
existence theorem (again, see [17]). that has been uniquely extended using Theorem 2. Since this may be 
done as long as the p-norm stays bounded, we will only need to prove the following lemma to establish 
Theorem 1. 

Lemma 1 Assume q > 7 + a/33 and conditions (I), (II) and (III) hold. Let f be a C 1 solution of (QP on 
[0, T] x R 3 x R 3 . Then, 



< C 



(F -/)(*) dv 

for all t € [0, T], where C is determined by F, A, fo, and T 



Although it is not explicitly stated in [17] , the proof presented there shows that <5 in Theorem 2 is bounded 
away from zero as long as || J(F — f)(t) dv\\ p is bounded. Using this observation, Theorem 1 follows from 
Lemma 1. 

Once Lemma 1 has been established, we will show increased decay of the charge density, p, under slightly 
modified assumptions. Instead of (III), we will assume the following conditions for some C > and all 
(>0,ieR 3 , and v <E JR 3 : 

(IV) A is C 1 with 

\A(t,x)\ <CR- 2 (x), (7) 



\d x A(t,x)\<CR- 3 (x), (8) 

and 

\V x -A(t,x)\=0. (9) 

Notice that condition (IV) does not involve spherical symmetry of A, and thus the result which follows 
requires no such assumption. 

Theorem 3 Assume conditions (I) and (II) hold, without ([5j) ; and condition (IV) holds. Let T > and f 
be the C 1 solution of (QP on [0,T) x M 3 x R 3 with 

|||(F-/)(t)|||<oo 

for every t £ [0,T). Then, we have 



I|p(*)lla < C Ptt 

for any t £ [0, T), where C p j depends upon 

sup \\p(r)\\ P - 

r£[0,t] 



In Section 1 we will prove Lemma 1, and thus Theorem 1. Then, Section 2 will contain the proof of 
Theorem 3. We will denote by "C" a generic constant which changes from line to line and may depend upon 
/o, A, F, or T, but not on t, x, or v. When it is necessary to refer to a generic constant which depends 
upon other variables, we will use variable subscripts to distinguish them. For example, we will use C Pjt quite 
frequently in Section 3, and it will always denote dependence upon ||p(t)|| p and t. When it is necessary to 
refer to a specific constant, we will use numeric superscripts to distinguish them. For example, will 
always refer to the same constant. 



Section 1 

Define the characteristics X(s, t, x, v) and V(s, t, x, v) by 



Then, we have 



^■(s,t,x,v) = V(s,t,x,v) 
f£(s, t, x,v) = - (E{s, X{s, t, x, v)) + A(s, X(s, t, x, «))) 
X(t, t, x, v) = x 
V(t, t, x, v) = v. 



(10) 



—f(s, X(s, t, x, v), V(s, t, x, v)) =d t f + V- V x f — (E + A)- V v f = 0. 
ds 



Therefore, / is constant along characteristics, and 



f(t, x, v) = /(0, X(0, t, x, v),V(0, t, x, v)) = f o (X(0, t, x, v), V(0, t, x, v)). 



(11) 



Thus, we find by (17) that / is nonnegative and sup^. / = ||/o||oo < oo. Unless necessary, we will omit 
writing the dependence of X(s) and V(s) on t, x, and v for the remainder of the paper. 

In order to bound the electric field, and thus the velocity support, we must use Lemma 3 and Theorem 
4 from |16j . In particular, we state the following lemma without proof. 

Lemma 2 Assuming conditions (I) and {III) hold, without (iSj), there exists k : [0, T] x M 3 — > [0, oo) such 
that 

\p(t,x)\ < C(k(t,x)i +k(t,x)i) 



and 



k(t,x)dx < C. 



Furthermore, notice that, due to the radial symmetry, we may write 



where the enclosed charge is given by 



= m M g|) £ (12) 
\xr \x\ 



m(t,r) := / p(t,y)dy. 

J\y\<r 



Lemma 2 will be exactly what is needed to bound the electric field, E. We have for any x G 



/ k(t,y)$dy < (f k(t,y)dy) (f dy 
J\v\<\x\ \J\y\<\*\ J \J\y\<\x\ , 



< / k(t,y)dy) (C\x\ 3 y 
\J\y\<\x\ J 

< C\x\i 



and 



3 

/ k(t,y)idy < I f k(t,y)dy) (C\x\^ 
Av\<\x\ \ J \y\<\*\ J 

< C\x\i. 



Now, we use the bounds on k to estimate the enclosed charge. 



\m(t,\x\)\ < [ \p(t,y)\dy 



Thus, 



So, 

and combining this with (III), 
where 



< C (k(t,y)T +k(t,y)*) dy 

J\y\<\x\ 

< C(\x\% + \x\i). 



\E{t,x)\ < C\m{t,\x\)\\x\- 2 

< C\x\- 2 (\x\i + \x\i) 

< C(\x\~i + \x\-*). 

\E(t,x)\<CQ(\x\), 
\E{t,x)+A(t,x)\ <Cg(\x\) 



r s , r < 1 
r~ 5 , r > 1. 



Next, we use methods from [7] to estimate the velocity support. Assume 

f(t, x, v) = /(0, X(0, i, x, v), V(0, t, x, «)) ^ 
so that, using the compact support in v of /o, we have 

|X(0)| = |V(0)|<C 
For any e > 0, define a := 2 + e, b :— -^zlj, and 



B:= (2 J™ G(\x\) a dx)i. 



Let 4i,t2 € [0,T] with t\ < £ 2 , and assume Xi > on [ti,t 2 ]. 
We claim 



/•X»(t 2 ) 

/ SOI) cte < fl(Xi(ta) - Xi(ti))v + 10. 
^i(ti) 

Si (a:) := 



To show (HHJl, let 



and 



Then, 



* 4 (t 3 ) 



Xi(tl) 



£7(|jc|) dx 



X 2 



Gi(x) dx + Q2{x) dx 



if |a:| < 1 

if \x\ > 1 

if \x\ < 1 

if \x\ > 1. 



x 4 (t 2 ) 



< 



(?i(x) da; 



Xi(tj) 



(£ 2 (x)) Q (^(tiJ-JCi 



-1 \JXi(h) J 

/ poo \ a 

< 10+ ( 2^ (£ 2 (x)J daO° (X^-X^))* 

< 10 + S(X 1 (t 1 )-X l (< 2 ))i 

This establishes (|13|) . Now, using this result and following [7j, we have 



\Xi(t 2 ) 2 - Mh) 2 \ = 



1-2 



2 / Xi{a)Xi{a) ds 
Jti 

< C f 2 G(\X(s)\)X t (s) ds 

< C ( 2 g(\Xi(a)\)Xi(a) ds 

Jti 

=: / G(\x\) dx 

JXi(h) 

< C^B(Xi(t 2 ) - Xi(ti))i + 10C<°> 

< C^B{{ sup \Xi(T)\)(t2-ti))i+10CW. 

re[*i,ta] 



Define 

W:= sup \Xi(s)\. 
se[o,T] 

Then, 

\X(t 2 ) 2 - *<(ii) 2 | < C^B(Wt 2 )* + lOCW. (14) 
Note that this holds if X; < on [ii,i2] 5 as well. Let us consider t G [0, T] and X(i) > 0. Define 

t:= inf-jr > : X*(s) > 0, Vs G [r,f]}. 

If t = 0, then by QX 



< Xf(0) + C^B(Wt)i + 10C* (0) 

< X?(0) + C(°>5(iyT)i + lOC^. 

If i > 0, then Xj(i) = and by l[I3|). we have 

< Xf(t) + C'( 0) B(Wt)* +10C (0) 

= c (0) s(i^r)^ + ioc* (0) 

< |Xf(0)| + C<°>.B(WT)i + 10C (0) . 
We may repeat this process for X,(i) < 0, so summing over i yields 

\X(t)\ 2 < |X(0)| 2 + 3C^B(WT)i + 30C (0) . 
Since the right-hand side is independent of t, we take the supremum and find, 

W 2 < ZC {0) BT^W^ +M 
where M := |X(0)| 2 + 30CW. If 3C^BT^W^ < M, then 

W < V2M =: d(T). 

If M < SC^BTIw*, then 

W < (6C {0) B)^T^ := C 2 (T). 

So, combining the inequalities 

W < max{Ci(T),C 2 (T)}. 

Thus, all characteristics V(s,t,x,v), along which / is non-zero, are bounded for any s G [O,? 1 ], including 
V(t, t, x, v) = v. Define 

Qt := sup{|w| : 3i e I 3 , t G [0,t]such that/(r,x,w) ^ 0}. 
Notice that Q is an increasing function of time, so that we may write for every s G [0, T] 

\V(s,t,x,v)\ <Q S < Qt- 
Since the momentum is bounded, bounds on position follow. Note that 

\X(s)-x\ < J\v{r)\dr 
< QtT. 



So, for |x| > 2Q T T, 



and 



\X( S )\>\x\-Q T T>\\x\ 



\X(s)\<\x\+Q T T<^\x\. 
Hence, we have control on X characteristics. 

Now, to bound the charge, we must first estimate the corresponding density. We use ([2]) to find 

±(g(s,X(s),V(s))) = -(E( S ,X( S )) + A(s,X(s))).V v F(V(s)). 
as 

Thus, 

rt 



g(t,x,v)=g(0,X(0),V(0))- / (E(t,X(t)) + A(t,X(t))) ■ V v F{V{r))d T 

Jo 

and by (III), @2J), and 

\g(t,x,v)\ < \ 9o (X(0,t,x,v),V(0,t,x,v))\ 

+ J*(\E(T,X(T,t,x,v))\ + \A(T,X(T,t,x,v))\) \V v F(V(T,t,x,v))\ dr 



< 



C\X(0,t,x,v)\- 2 + C f*\X( T ,t,x,v)\~ 2 HV^Hoo dr 
+C f*\m{T,X(T,t,x,v))\ |A:(r,t,a:,t;)|- 2 ||V t ,f , || o dr 
< C\x\- 2 (l + J*\m(T,X(T,t,x,v))\ dr) 

for \x\ > 2Q T T. 

To show that g decays like r~ 2 , we must bound the enclosed charge m. For |x| < 2QtT, 



\m(t,\x\)\ < I \p(t,y)\dy 

\y\<2Q T T 



< (2Q T T) 3 [ ||F-/o| 

J\v\<Qt 



dv 



C 



Now, let r = \x\ and define M.{t) := sup r \m(t,r)\. We know by ©, (fT5|) . and the Divergence Theorem, 



\m(t,r)\ < \m(0,r)\ + 
= |m(0,r)| + 



d s (g(s,y,v)) dv dy ds 



'0 J\y\<r J\v\<Q T 



' l + / 'M(r) 

\y? 



dv dS y ds 



v ■ V y g{s, y, v) dv dy ds 

J\y\<r J\v\<Q T 

< |m(0,r)|+/ / / H \g(s,\y\,Kyv)\ dv dS v ds 

JO J \y\=r J \v\<Q T 

< M(0) + C f f f \v 

JO J\y\=r J\v\<Q T 

= M(0) + C f Q$.(l+ [ M(r) dr) ds 
Jo Jo 

< C + C f M(t) dr 

Jo 

for r > 2Q T T. Then, since m(t, r) < C for r < 2Q T T, we find, for all r, 

m(t,r) <C + C [ M(t) dr 
Jo 



and consequently 



M{t) <C + C M(t) dr. 
Jo 



By Gronwall's Inequality, 

M{t) < (M(0) + CT)(1 + Ctexp(Ct)) < (M(0) + CT)(1 + CTcxp(CT)) < C. 
Notice, then, that this bound and (jT5J) imply 

C 

\g{t,x,v)\ < 

for |i| > 2Q T T. 

Proceeding in the standard way, we next estimate the gradient of the electric field. 



(19) 



d a ,E k = 



OXi ylx] 6 J\y\<\ x \ J 



dx 



XI 



3 dxj 



Then, 



and, letting r = \x\ 



d ( x k \ 






dxi \\x\y 




\x\* \x\ 5 



< 



c 



—m{t,r) 
ox t 



(*>»•) (7)| <C\p{t,r)\i 



Thus, 



mm t m C C|p(t,r)| r 2 



< 



C|p(t,r)|. 



Since the best we can do a priori is \p(t, r)| < Cr 2 , we can only conclude 

\d Xi E k \ < C\x\- 2 . (20) 

We begin to estimate the spatial decay of p by first estimating the large |x| behavior of the field integral 
obtained by integrating (| 1 7|) in v. Assume |x| > 8QtT and f(t,x,v) is nonzero. Define 



£(t,x) = E(t,x) + A(t,x). 



Then, 



£(s,X(s))-V v F(V(s))dv 



\v\<Qi 



< 



£(s,X(s)) ■ (V v F(V(s)) - V v F(v)) dv 



v\<Qt 



(£(s, X(s)) - £(s, x+(s- t)v)) ■ V v F(v) dv 



v\<Qi 



V„ • (F(v) 8(s, x+(s- t)v)) dv 



v\<Q T 



F(v)V v ■ (8(s,x+ (s - t)v)) dv 

\v\<Qt 

I + 11 + III + IV. 



By the Mean Value Theorem, (12]). ifTS]). and CE 



I < I {v1 < Qt \8(s,X(s))\ IIV^IU \V(s)-v\dv 

< ij^^CClX^I^HV^IU \f s E(T,X(T)) dr\) dv 

< CQ 3 T \x\- 2 (t-s)\x\- 2 

< C\x\~ 4 . 



(21) 



To estimate 77, we again use the Mean Value Theorem and (f20|) . so that there is £ x between X(s) and 
x + (s — t)v with 



Thus, wc find 



\8 i (s,X(s))-8 i (s,x+(s-t)v)\ = \V x £i(s,£ x )-(X(s)-x-(s-t)v) 

< C\Z x \- 2 \X(s)-x + (t-s)v\. 



II < C I \£ x \~ 2 \X(s)-(x + (s-t)v)\ ||V„F||oo dv 

l\v\<Q T 
J\v\<Q T 

< CT 2 Q 3 T \x\- 2 \£ x \- 2 , 



\V(t) -v\d,T] dv 



and by (flip)) . 



|&| > \X(s)\-\X(s)-(x + (s-t)v)\ 



> 5 N 



(7(r) - v) dr 



> -|z|-2TQ T 



> jM 



Therefore, 

By the Divergence Theorem, we find 
Then, evaluating 7V yields 



//< C|a;|- 4 . 

m = o. 



Since 



we have 



IV < / \\F\\ 00 \V x -£(s,x + {8-t)v)\\s-t\dv 

J\v\<Q T 

< CT \p(s,x+(s-t)v)+V x -A(s,x + (s-t)v)\dv 

J\v\<Q T 

< C I \p(s,x+(s-t)v)\ dv + CR- i (x+{s-t)v). 

J\v\<Q T 



\x + (s-t)v\ > \x\-TQ T >-\x\, 



IV <C I \p{s,x + (s-t)v)\ dv + C\x\^. 

l\v\<Q T 



Finally, collecting (J2IJ, (EH), flSSJ), ©, we have for |x| > 8Q T T, 



£(s,X(s)) • V„F(V(*)) <fo 



<C|a;r 4 + C/ |p(«,x+(«-t)«)| dv. 
\v\<Qi 



\v\<Q T 

Now, we can bound ||/?(f)||4. Using the bound on the velocity support, we have 



(22) 
(23) 



(24) 



(25) 



l|p(*)ll 



sup 

X 

4tt 



\v\<Q t 



(F(v)-f(t,x,v)) dv 



Thus, for \x\ < 8Q T T, we have 



C. 



iP(*,x)i<iiP(t)iioo <ci*r 4 . 



Now, recall ffTTj) 



g(t,x,v) = g(0,X[0,t,x,v),V(0,t,x,v)) - / £(s,X(s,t,x,v)) -V v F(V(s,t,x,v)) ds 



We use this equation and (j2"5)l so that for \x\ > 8QtT, 



\p(t,x)\ 



g(t, x, v) dv 



v\<Q T 



< I ( \g o (X(0),V(0))\ 

! \v\<Qi 



S(s,X(s))-V v F(V(s)) ds 



dv 



< CQ$\x\ 



< C\x\~ 4 + C 



[ [CQ 3 T \x\- 4 + C I 

JO \ A v \<Qt 



\p(s, x + (s ~ t)v)\ dv j ds 
\p(s, x + (s — t)v)\ dv ds. 



o J\v\<Q T 



So, 



Define 



Then, for all x £ 



\x\ 4 \p(t,x)\<C + C / \x\ A \p(s,x + (s-t)v)\ dv ds. 

JO J\v\<Q T 

V(t) :=sup(|x| 4 \p(t,x)\). 



V{t) < C + C f f T(s) ds dv 
J\v\<q t Jo 
t 



< C + C f T(s) ds 
Jo 

and using the Gronwall Inequality, we find 



V{t) < C 

and thus \p(t, x)\ < C\x\~ 4 for all x G M 3 . Finally, we may apply Theorem 2 with any q > 7 + \/33, and 
since ||p(£)||oo is bounded, we conclude that \p(t,x)\ < CR~ 4 (x). Since this estimate is independent of T, 
we find 

sup \\p(t)\\ p < C. 
te[o,T] 

This remains true for any T > 0, so the proof of Lemma 1, and thus Theorem 1, is complete. 



Section 2 

As in Section 1, we will show that the charge density decays at a faster rate than previously known. In 
the work that follows, we will use the framework of the previous sections and assume conditions (J) and 
(//) from Section 1, without the spherical symmetry of /o- However, we will not take condition (///) as an 
assumption, and instead, assume condition (IV) holds for some C > and all t > 0, x € M 3 , and v £ M 3 . 
Since we have made a change in the assumptions, we may not use results from the previous sections, unless 
otherwise stated. Thus, this section can be viewed independently from the others, as we will rely more 
on results shown previously in [17] . Recall from the introduction that we will use C Pj t to denote a generic 



constant which depends upon ||/o(i)|| p and t. 

To begin, we apply Theorems 1 and 2 of [17] . finding a unique / G C 1 ([0,T] x R 3 x R 3 ) which satisfies 
(U} and assume 

|||(F-/)(t)|||<oo (26) 
for all t G [0,T). To first bound the velocity support, we write, as before, 

g(t,x,v) := F(v) - f(t,x,v) 

and 

£(t,x) := J5(t,ar) + A(t,a;). 
Using Lemma 1 of [17] . we find for all i£l 3 and t G [0,T), 

' \S(t,x)\ dr<C sup ||p(r)|| p . 



r£[0,t] 



\V(s)-V(0)\ < / |f(r,X(T,t, dr 



Then, for any s G [0,t), 



Assuming / ^ 0, we know f(t,x,v) — fo(X(0), V(0)), and by (II) it must follow that 

\V(0,t,x,v)\<C. 

Thus, |V(s)| < C PiS for any s G [0,i). Following previous notation let us write 

Qt := sup{\v\ :3x eR 3 ,T £ [0,t) such that f(r,x,v) ^ 0} 

so that for all s G [0,i) 

< Qt- (27) 

This bound on the velocity establishes some of the relations shown in previous sections. Most importantly, 
(fT5| and ^ must hold for |x| > 2Q t t. 

Next, we denote the ^-derivatives of characteristics by 

Ay := -r^- (s,t,x,v) 



and 



Then, using the characteristic equations of ([28 



Thus, 



dV 

——(s,t,x,v). 



ff(s)=B(.s), 
':(s) = V x £(s,X(s)) A(«), 
A| s=t = 0, 
B| s=t - I. 



|A(*)| < / |B(r)| dr 



(28) 



and 



< 1+ f \W x £(t,X(t))\ |A(r)| dr 

J s 



which leads to 



\A(s)\ + \M(s)\ <l + J (|B(r)| + |V x £(r,X(r))| |A(r)|) dr. 

Again, using Lemma 1 from [17j . we find 

\V x £(t,X(t))\ <C p , t R- 3 (X(t)). 

So, define 

"H(s) := sup (|A(s, i, a;, v)\ + |B(s, i, 

{(x,v):\x\>2Q t t} 

Then, for |x| > 2Q t t, we use (HHJ) to find 

R- 3 (X(t)) < R- 3 (Q t t), 

and thus, from (|29[) 

H(s) < l + max{l, C p , T i?^ 3 (Q t t)} W(r) dr. 
Using the Gronwall Inequality, we find 

< C p<t 

for s G [0,t], and v-derivatives of both characteristics are bounded for \x\ > 2Q t t. 
Summarizing ()28|) , we may write 

^(s)=V x £(s,X(s))A(s); A\ s=t =0; ^jUt=L 



So, 



Again using 



Thus, 
where 

Now, 



A(a) = (s-t)I + J J V X £(X,X(X))A(X) dXdr 
=: (s — i)I + 7i(s, t, x, v). 



— (s) = V x £( S ,X( s ))A(s) 

= \7 x £(s 7 X(s)){{s - t)l + 71 (s, t, x, »)) 
=: (s - t)V s £(s, X(s)) + 72 (s, i, x, «). 

M(s) = l + l3 (s,t,x,v) 
j 3 (s,t,x,v) = - f ((t - t)V x £(T,X(r)) +~f 2 (T,t,x,v)) dr. 



— (detB) = det E 21 B 22 B 23 ] + dct ( B 21 B 22 B 23 ] + det 

> ^5 TOJ TK 

*31 B 32 Ji»33 
I + 11 + III. 



'11 


"12 


JJJ 13 


•21 


B 22 


B23 


'31 


B32 


B33 



Estimating the first term yields 



Similarly, 



and 



So, we find 



and since 



(s-*)f& + (72)n ( S -«)fg; + (7 2 )i2 (a-t)^ + (72)13 

det ( (73)21 1 + (73)22 (73)23 

(73)31 (73)32 1 + (73)33 

go 

h (72)11 +<r 1 (s,t, a;, v). 

oxi 



d£ 

II =■ (s - t)- h (72)22 + cr 2 (s, t, x, v) 

0x2 

d£ 3 

III =: (s - t)- 1- (72)33 + (r 3 (s,t,x,v). 

dx 3 



d 3 
— (det B) = (s - t)V x ■ £(s, x)\ x=x{s} + ^((72)^ + <Tj 

3=1 



detB = 1 - 4tt / (r- f)p(r,X(r)) dr - / ^((72)^ + ° 3 ) dr 



Let 

r t 3 



e{s,t,x,v):=4:TT (t -t)p(r,X(j)) dr + ^(Mjj + dj) dr. 

J S J S 

Then, for |e| < 1, 



1 

l-e 



= i+c+Er^ n (33) 

=: 1 + 4tt f*(r-t)p(T,X(T)) dr + r)(s,t,x,v). 

Now that the determinant has been written in a nicer form, we estimate the remaining terms. Let 
\x\ > 2Q t t. Using (fT5|) . (|26|) . (|30| . ff3T~|) . and (|32|). we find the following bounds on the error terms for any 
i,j = 1,2,3: 

|(7i)ij| < C p , t T 2 i?- 3 (x) sup |A(s,t,z,«)| 

se[o,t] 

< C p , t R- 3 (x), 

1(72)^1 < C p , t i?- 6 (x), 

and 



1(73)^1 < c^r^^ + c^i?- 6 ^) 

< C p , t ir 3 (x). 



Then, using the estimates of the 7 terms, for any k — 1,2,3, 

k fc | < C(|V :c f| + |72|)(2|73| + 2|73| 2 ) + 2(|V :!; f| + | 72 |)(|73| + 2|73| 2 ) 



< [C p . t TR- 3 (x) + C p , t R- 6 (x)}[2C Pit R- 3 (x) + 2C p , t R- 6 (x)} 



+ 2[C p . t TR- 3 {x) + C p , t R- 6 }[C p , t R~ 3 {x) + 2C p , t R~ 6 (x)} 

< C p , t R- 6 {x). 
Now, we may bound e. 

|e| < AttT 2 sup \\p(T)\\ p R- p (x)+ 3T sup (| 72 (r)| + |a(r)|) 



r£[0,t] 

< C Pit R- p (x) + C p>t R-*(x) 

< C p , t R~ p (x) 
=: C (1) R- p (x) 



r€J,t 



1 

< 2 



for |x| > (2C^ 1 ^)p . Thus E^L 2 e ™ converges and rj is well-defined for large enough values of Finally, we 
have, for \x\ > max{2Q t t, (2C«) * }, 



\v\ 



< 3C p , t Ti?- 6 (a;) + ((C«) 2 ir 2 P(x) + (C^^ir^x) + ...) 

< C p , t R~ 6 (x). 
Using ((271), we find for t € [0,T), 

||p(t)||oc = sup f(F(v)-f(t,x,v)) dv 

x J 

< -|q 3 \\f -/oiu 
— ^p,<- 

Thus, if |x| < .D for some D > 0, we have 

|p(t,i)| < \\ P (t)\\ooR 6 (D)R- e (x) 
< C p , t R~ 6 {x). 

Now, denote := max{2Q t i, (2CW) * , 27V}, and let |x| > . Then, by JS) we have 



(34) 



(35) 



\X(0,t,x,v)\ > -\x\ > N. 



Using (V), ([IT]). ([55)1. and fl5J} to estimate p(i,x) yields 
/(*, x, v) dv = I f (X(0, t, x, v), V(0, t, x 7 v)) dv 

F(V(0,t,x,v)) dv 

F(V(0,t,x,»))(detB(0)) 5 ^ 5 r dv 



, * dV 
F(V(0, t, x, v))(l + 4?r / (r - i)p(r, X(t, t, x, v))dr + r)(0, t, x, w))(det ^—(0, t, x, v)) dv 



dv 



dV 

F(w) dw+ I F(V(0,t,x,v))ri(0,t,x,v)(det—(0,t,x,v)) dv 



f /"* dV 

+4tt J F(V{0,t,x,v)) J (r -t)p(T,X(T,t,x,v)) dr (det — (0, t, x, v)) dv. 

Since p(t,x) = J (F(v) — f(t,x,v)) dv, we have 

p{t, x) = Am j F(V(0, t, x, v)) J*(t - r)p(r, X(r, t, x, v)) dr (det §£(0, t, x, «)) ^ 

(36) 

- / F(V(0, t, x, v)) V (0, t, x, v)(det §£(0, t, i, w)) dv. 

Now, let 

tf(t) := ||p(t)|| 6 = sup(|p(t, ar)j^ e (^)). 

In order to make the change of variables w — V(0,t,x,v) in the v-integral, we must first show that the 
mapping v — > V(0, t, x, v) is bijective. For the moment, we will take this for granted, and continue with the 
estimate of p, delaying the proof of this fact until the very end. By ([34| . ([36)l. and the work of Section 3.6, 
we have for \x\ > C^ 2 \ 

\p(t,x)\ < 4.TrfF(V(0,t,x,v))(f*(t-T)V(T)R- 6 (X(T,t,x,v)) dr)(det%(0,t,x,v)) dv 
+C p ^tR' 6 (x) J F(w) dw 

< CR~ 6 (x)fF(V(0,t,x,v))(det^(0,t,x,v))dvf<(t-T)y(T)dT (37) 

+C p jR- 6 (x) 

< R- 6 (x)(C / *(*-t)*(t) dT + C p , t ). 

So, applying ()35l) with D = and combining with ([57| . we have for all x, 

\p(t, x)\R 6 (x) < C p ,t +C [ (t - r)*(r)dr 

Jo 

and since the right side is independent of |x|, 

V(t)<C Ptt +C f (t-T)*(r)dT. 
Jo 

By the Gronwall Inequality, 

*(*) < C*,t- 

Therefore, for every f £ [0,T), 

l|p(*)lla < C Pl( 

and the proof is complete. 



Change of Variables 



In order to justify the change of variables used in line (1371) of the previous section, we must first demonstrate 
that the mapping v — > V(0, t, x, v) is bijective. This is done below. 

From (|30|) . we know there is C^fl > such that 



\W£{t,x)\<cf}\x 



for every x G M 3 , t G [0,T]. Also, using ([32]), we know there is > such that for all s,i£ [0,T], 
x,v G M 3 , with |x| > 2TQ(T), 



— (s,t,x,v) 

(5) 



< c, 



(4) 



P,* • 



Finally, using Lemma 1 of [17], we know there is Cp t . > such that for all a; £ M 3 and t £ [0,T], 

|£(r,x)| dr < C$. 

For any D > 0, define C D := max{8 (.D + C$T\ ,2TQ(T),4(6C$C$T)i} and B(0,£>) := {v G R 3 : 
< -D}. Injectivity on B(0, 1?) can now be shown in the following lemma. 

Lemma 3 For any D > 0, \x\ > Cd, and t G [0,T], the mapping v — > V(0,t,x, v) is injective on B(0,D). 

Proof: Let D > 0, |x| > Cd and t G [0,T] be given. Then, let wi,u 2 G B(0,D) be given with 

V(0,t,x,»i) = V(0,t,a;,«a). 

We have 

vi+ I £(t, X(r,t,x,vi)) dr — v 2 + / £ (r, X(t, f, x, U2)) dr. 
Jo Jo 

So, using the Mean Value Theorem, for any i = 1,2,3 there is between X(t, t, x, v±) and X(t, t, x, U2) and 
between «i and «2 such that 



and 

Then, since 



£i(T,X(T,t,x,Vi)) - £i(T,X(T,t,x,v 2 )) = V x £i{T,C x ) ■ (X(r,t,x,vi) - X(T,t,x,v 2 )) 
Xi(r,t,x,V!) - X l (T,t,x,v 2 ) = V v Xi(T,t,x,C v ) ■ (v 2 - ui). 





1 




> 




x|- 




2 




1 




> 




x|- 




2 




1 




> 




x\- 




2 




1 




> 




x\- 




2 




1 




> 




x\. 




4 



|£| > |A:(r,t,a:,«i)|-|X(r,t,a;,t;i)-X(r,t,a;,«2)| 
"I _ / \V(ht, x, vi) - V(t,t, x, vz)\ dc 

ct 



v 1 -v 2 \ + 2J J ||-E(A)||oo dX dt 



and 



for any i = 1,2, 3, we find 



\Xi(T,t,X,Vt) - Xi(T,t,X,V 2 )\ < C^J \v 2 - Vi\, 



\Ei(T, X{t, t, x, Vl )) - £i{r, X{t, t, x, v 2 ))\ < |V k ^(t, £)\ \X(t, t, x, Vl ) - X(t, t, x, v 2 )\ 



< 



^supC^|^|- 3 j (>/3C$ \ Vl -v 2 \) 
< ^}c^(]\x\]\ Vl -v 2 \. 



Therefore, we have 



\v\-v 2 \ < / \£(T,X(T,t,x,V!))- £(T,X(T,t,x,v 2 ))\ dr 
Jo 

•ah 



rt /i \ -3 



□ 



< 3<7<>i-t»| J \j\x\ ) dr. 

< {ZCf}C$T{\\ x \)-*)\ Vl -v 2 \ 
1, 

< -j\vi-v 2 \. 

Thus, \v\ — v 2 \ =0, which implies v\ — v 2 , and injectivity follows. 



Now, let t e [0,T] and x e K 3 be given. Define 

S := {w : F(w) ^ 0} 

and 

V-^S) :={v:F(V(0,t,x,v))^0}. 
Using the compact support of F, we conclude that v e V~ 1 (S) implies 

\v\<W+ f 'WEWWn dr < c p>t . 
Jo 

Therefore, there is a D > 0, such that U _1 (5) C B(0,D). Thus, for |x| > C D , and t E [0,T], the mapping 
v — > U(0,t,x,v) is injective on V~ 1 (S) and bijective from V~ 1 (S) to 5. Finally, 

J F(V{0,t,x,v))det (jj^J dv = J F(V{0, t, x, v)) det (j^j dv 

= / F(w) dw 
Js 

= J F(w) dw. 

Thus, the change of variables is valid and the justification is complete. 
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